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Objective The risks of cigarette smoking concerning higher systemic disease mortality are lessened by smoking cessation. 

Methods Microarray analysis compared the expression profiles of smokers who were successful and not successful at smoking ces- 
sation, with the goal of identifying genes that might serve as potential biomarkers or that might be valuable in elucidating distinct bi- 
ological mechanisms. The mRNAs were isolated and compared from peripheral leukocytes of six smokers who were successful in 
cessation and six smokers who failed in smoking cessation. 

Results Two hundred ninety nine genes displayed significantly different expression; 196 genes were up-regulated and 103 genes 
were down-regulated in the success group compared to the failure group. Twenty four of these genes were identified with biological 
processes including immunity, cytoskeleton and cell growth/cycle. Real-time PCR confirmed the differential gene expression. The 
mRNA levels of HEPACAM family member 2 (HEPACAM2) and tropomodulin 1 (TMOD1) were significantly more expressed in 
the success group, while the mRNA ubiquitin specific peptides 18 (USP18) were significantly less expressed in the success group 
compared to the failure group. 

Conclusion The results suggest that smoking cessation can modulate cell adhesion and immune response by regulating expression 
levels of genes, especially HEPACAM2, TMOD1 and USP18, which have an important relationship with smoking cessation. 
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INTRODUCTION 

Cigarette smoke profoundly affects the development and 
progression of many diseases, including cancers, cardiovascu- 
lar, chronic obstructive pulmonary and gastrointestinal diseas- 
es. 1,2 Studies addressing the various adverse effects of smok- 
ing, including carcinogenic processes, vasomotor dysfunction, 
impaired endothelial-dependent vasodilatation, and lipid mod- 
ifications, 3 have established that smoking cessation is crucial- 
ly important in lessening and preventing smoking associated 
preventable diseases and promoting improved public health. 
Recent studies addressing the effectiveness of smoking cessa- 
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tion have shown that quitting smoking improves the overall 
life expectancy and lowers the risk of lung cancer, heart attack 
and stroke. 4 

Cigarette smoke is the principal environmental risk factor. 
Cigarette smoke negatively influences the immune system, es- 
pecially interferon induction, apoptosis and xenobiotic me- 
tabolism-related gene developments, and is important in neu- 
trophil-derived oxidant injury. 5,6 Stopping smoking is associated 
with the elevation in collagenase mRNA expression in human 
neutrophils, which are crucial in host defense against micro- 
organisms as primary phagocytic cells, and in acute phases 
of inflammatory reactions. 7 Although several genes have been 
linked to the incidence of smoking-related diseases, 8 little is 
known of gene expression profiling as an alternative approach 
to discover positional candidates of benefit to address dam- 
age attributed to smoking cessation. 

In the present study, microarray analysis clarified the gene 
expression profiles of peripheral leukocytes from smokers who 
succeeded or failed in smoking cessation. The differential ex- 
pression of the selected genes in smokers who were successful 



290 Copyright © 2014 Korean Neuropsychiatry Association 



SJ Kim et al. 



or unsuccessful in quitting was confirmed using real-time poly- 
merase chain reaction (PCR). 

METHODS 
Participants 

Twelve male subjects (mean age 30.4±6.8 years) with nicotine 
dependence defined according to the Diagnostic and Statistical 
Manual of Mental Disorders, Fourth Edition, Text Revision 
(DSM-IV-TR) criteria were recruited at Bucheon St. Mary's 
Hospital, The Catholic University of Korea, by advertisements 
in the hospital newspaper and on websites. Participants had 
no mental or physical afflictions, aside from their nicotine 
dependence. They were not using illegal drugs and any pre- 
scribed medications. The participants were instructed not to 
drink more than five standard alcoholic drinks per week and 
one cup of coffee per day, and to maintain their eating habits 
during the smoking cessation to minimize the influence of al- 
cohol, caffeine and change of diet habits on mRNA expres- 
sion levels. 

Smoking cessation 

During the current study, all subjects agreed to quit smok- 
ing without taking medication or nicotine replacement thera- 
py. Smoking abstinence was assessed from self-report. Urine 
cotinine levels were tested by a urinary cotinine immunoassay 
dipstick (Nicometer; Jant Pharmaceutical, Encino, CA, USA) 
prior to smoking abstinence and one and two months after 
smoking abstinence. Success in smoking cessation was validat- 
ed by a urine cotinine level <30 ng/mL. All subjects signed an 
informed consent after receiving a complete description of 
the study. The study was approved by the Ethics Committee 
of Bucheon St. Mary's Hospital. The procedures were in ac- 
cordance with the Helsinki Declaration of 1975, as revised in 
1983. Six of the 12 smokers failed to abstain from smoking 
during the study (failure group), as revealed by the results of 
urinary tests. The remaining six smokers were successful in 
smoking abstinence (success group). 

RNA preparation 

After two months, fasting peripheral blood of each subject 
was drawn into an EDTA-coated tube treated with 10 uL/mL 
of Pefabloc SC (Roche Diagnostics GmbH, Mannheim, Ger- 
many). One milliliter of each suspension was mixed with 2 
mL of TRIzol LS reagent (Invitrogen, Carlsbad, CA, USA) and 
total RNA was isolated and purified using the RNeasy mini 
kit (Qiagen, Valencia, CA, USA). The amount of RNA was de- 
termined by measuring the optical density at 260 nm (OD 2 6o) 
with a Nanodrop apparatus (Thermo Scientific, Wilmington, 
DE, USA) and the quality of the RNA was measured using a 



model 2100 bioanalyzer (Agilent Technologies, Santa Clara, 
CA, USA). 

Expression profiling 

Two hundred nanograms of each total RNA was used to 
make Cy3-labeled cRNA using Low Input Quick Amp La- 
beling Kits (Agilent Technologies). The labeled cRNAs were 
applied to Custom Gene Expression Microarray, 4x44k (Agi- 
lent Technologies) at 65°C for 17 hours. After hybridization, 
the microarrays were washed according to the manufactur- 
er's recommended protocol. Microarray images were acquired 
using GenePix 4000B (Molecular Devices, Sunnyvale, CA, 
USA) and the signal intensity of each probe on each microar- 
ray was determined using Feature Expression Software (Agi- 
lent Technologies). Interarray variation was normalized by 
quantile normalization. After excluding weak signals (<10 2 3 ), 
the remaining 30,398 features were analyzed. 

Microarray data analysis 

To reveal differentially expressed genes between the smok- 
ing cessation success and failure groups, we calculated shrunk- 
en centroids of each group using the pamr package of R. 9 The 
maximum number of features showing no misclassification 
error was selected. 

Real-time PCR 

We conducted a RT-PCR to confirm the differential gene 
expression that selected HEPACAM2, TMOD1, USP18 and 
IF 144 between the success group (n=6) and the failure group 
(n=6). All oligonucleotide primers were synthesized by Bi- 
oneer (Daejeon, Korea) using the following primer sequenc- 
es: (3-actin-forward primer, TCA CCC ACA CTG TGC CCA 
TCT ACG A; (3-actin-reverse primer, GGA TGC CAC AGG 
ATT CCA TAC CCA; HEPACAM2-forward primer, CTG 
GGG CTG TGG AGT ATG TG; HEPACAM2-reverse prim- 
er, CAG GCA GCT GTA ATT CCC AA; TMODl-forward 
primer, TGA AGC CAC TGA ACT CTG CC; TMOD1 -re- 
verse primer, TTG TTG CCT GGG ATT GGT AA; USP18- 
forward primer, GTC CAT CCT GGC TGA GTC CT; USP18- 
reverse primer, ACC AGG CCA TGA GGG TAG TC; IFI44- 
forward primer, GCA TGT AAC GCA TCA GGC TT; and 
IFI44-reverse primer, CTC CTT TCT CAC TCA GCC CC. 
Real-time quantitative RT-PCR was performed in a 25 uL to- 
tal volume containing 20 ng total RNA, 7.5 uL RNase Free dis- 
tilled water, 1 uL forward and reverse primers (10 (iM), 1 uL 
PrimeScript 1-step Enzyme Mix 2 and 12.5 uL 2X One Step 
SYBR® RT-PCR Buffer 4 (Takara, Kyoto, Japan). The reaction 
conditions were set at 42°C for 5 min and 95°C for 10 sec, fol- 
lowed by 40 cycles of 95°C for 5 sec and 60°C for 30 sec. Real- 
time PCR was carried out using a thermal cycler CFX96 (Bio- 
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Rad, Hercules, CA) and triplicate 25 uL reactions were per- 
formed in Hard-Shell® PCR Plates 96- Well WHT/WHT (Bio- 
Rad). 

Comparative calculation 

The results of real-time PCR data are shown as Ct values, 
which are defined as the number of cycles required for the 
fluorescent signal to cross a threshold and which are inverse- 
ly proportional to the amount of target cDNA in the sample. 
The ACt value was calculated by subtracting the Ct of control 
(P-actin) from the Ct of target gene (HEPACAM2, TMOD1, 
USP18 and IFI44); A ACt represented the different expression 
of each sample, according to mean expression of the success 
group, as calculated by the formula AACt=ACt-mean ACt of 
the success group. The N-fold differential expression in a 
specific gene of a sample was expressed as 2 A (-AACt). 

Statistical analysis 

The comparison of the N-fold gene expression between 
the two groups, was performed by a two-tailed Students' t- 
test. The analysis was done using the Statistical Package for 
the Social Sciences (SPSS) version 18.0 for Windows. We con- 
sidered p<0.05 to be statistically significant. 

RESULTS 

Expression profiles 

We identified 299 features with no misclassification error 
using the pamr module of R (9). Hierarchical clustering based 
on the 299 features clearly divided the samples into two clus- 
ters; success vs. failure groups (Figure 1). However, individual 
differences hindered the finding of any common time-course 
change in each group. Of these 299 genes, 196 were up-regu- 



lated and 103 were down-regulated in the success group, com- 
pared to that of the failure group, and 24 genes were classified 
into biologically relevant categories according to their func- 
tions (Table 1). 

Real-time PCR 

To confirm the expression levels obtained by the microar- 
rays, relative quantitative real-time RT-PCR was performed 
for four selected genes: HEPACAM2, TMOD1, USP18 and 
IFI44 (Figure 2). Significantly higher expressions of HEPA- 
CAM2 (p=0.024) and TMOD1 (p<0.001) were evident, as was 
the reduced expression of USP18 (p=0.033) in the success 
group compared to the failure group. IFI44 expression was re- 
duced in the success group, but the decrease was not significant. 

DISCUSSION 

The present study used the microarray technique to explore 
the expression levels of genes as a way of assessing the effec- 
tiveness of smoking cessation in the peripheral leukocytes of 
the success group and the failure group. Expression of 299 
genes differed between the two groups, and 24 genes were 
classified into several functional groups. 

The first group included the genes involved in immunity and 
anti-immunity. APOBEC3A, CD55 were up-regulated and 
ITGB7, KIR2DS2, TNFRSF17, IFIT1, IFI44, IFI6, USP18, and 
HERC5 were down-regulated in the success group. The bio- 
logical effects of cigarette smoking cessation on the immune 
system are unclear. Smoking stimulates the oxidative burst, 
which has a potential role in disease progression by suppres- 
sion of neutrophil responses induced by subgingival plaque/ 
bacterial products and immune complexes disease sites. 10 ' 11 It 
is also involved in synthesis and secretion of immunoregulato- 




Figure 1. Hierarchical clustering of 299 
differentially expressed genes in chronic 
smokers. 
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Table 1 . Gene expression patterns of selective genes in the success group and the failure group 


Category 


Symbol 


Gene name 


Median expressior 


LOf 


Success (S) 


Failure (F) 


Ratio of S/F 


Immunity/ 


APU-DEC3A 


Apolipoprotein B mRNA editing enzyme, 


3.76 


2.84 


1.32 


Anti-immunity 




catalytic polypeptide-like 3A 












CD55 molecule, decay accelerating factor 
for complement 


3.39 


3.15 


l.Uo 






Integrin, beta 7 


4.05 


4.25 


0.95 






Killer cell immunoglobulin-like receptor, 
two domains, short cytoplasmic tail, 2 


2.55 


2.80 


0.91 




TNFRSF17 


Tumor necrosis factor receptor 
superfamily, member 17 


2.35 


2.68 


0.88 




MT1 


Interferon- induced protein with 
tetratricopeptide repeats 1 


3.82 


4.71 


0.81 




USP18 


Ubiquitin specific peptidase 18 


2.48 


3.18 


0.78 




HERC5 


Hect domain and RLD 5 


2.64 


3.61 


0.73 


Cytoskeleton 


TM0D1 


Tropomodulin 1 


3.06 


2.60 


1.18 




EPB49 


Erythrocyte membrane protein band 4.9 


3.32 


3.05 


1.09 




ANK1 


Ankyrin 1, erythrocytic 


3.34 


3.21 


1.04 




EPB41 


Erythrocyte membrane protein band 4. 1 


2.64 


2.53 


1.04 


Cell growth/Cycle 


MAPK3 


Mitogen-activated protein kinase 3 


2.80 


2.50 


1.12 




CCNG2 


Cyclin G2 


2.34 


2.10 


1.11 




ANXA3 


Annexin A3 


3.23 


3.08 


1.05 




IFI44 


Interferon- induced protein 44 


3.46 


4.32 


0.80 


Unclassified 












Antioxidant 


GSTT1 


Glutathione S-transferase theta 1 


1.99 


1.66 


1.20 


Calcium homeostasis 


NUCB1 


Nucleobindin 1 


2.69 


2.54 


1.06 


Cancer-related 


MXRA7 


Matrix-remodelling associated 7 


3.46 


3.66 


0.95 


Chaperone 


AHSP 


Alpha hemoglobin stabilizing protein 


4.13 


3.99 


1.03 


Development 


GRHL1 


Grainyhead-like 1 


2.38 


2.18 


1.09 


Digestion 


MGAM 


Maltase-glucoamylase 


2.88 


2.63 


1.09 


Protects against 


PINK1 


PTEN induced putative kinase 1 


3.80 


3.44 


1.10 


mitochondrial dysfunction 












Cell adhesion 


HEPACAM2 


HEPACAM family member 2 


2.40 


2.02 


1.18 



CD55: Complement decay-accelerating factor molecule, HERC5: Hect domain and RCCl-like domain-containing 
phatase and tensin homolog-induced putative kinase 1, HEPACAM2: Herpa cell adhesion molecule family member 



protein 5, PINK1: Phos- 
2 



ry cytokines and inflammatory mediators, and may reduce 
neutrophil chemotaxis that indicates alterations in the mor- 
phology and adhesion property of neurophils, which reduces 
functional capacity. 12 APOBEC3A prevents virus infections by 
providing intracellular resistance to retrovirus replication and 
transposition of endogenous retroelements. 13 The CD55 reg- 
ulatory protein of the complement cascade inhibits the net in- 
crease in complement-mediated cytotoxicity by the alternative 
pathway. 14 ITGB plays a critical role in the mucosal immunity 
and is linked to inflammatory bowel disease. 15 Additionally, 
KIR2DS2 influences natural killer cell activating, which pro- 
motes cell lysis and absence of KIR2DS2, and is associated 



with resistance to inflammatory bowel disease. 16 TNFRSF17 
might contribute to the progression of multiple myeloma 17 and 
IFIT1 is a potential diagnostic marker for myelodysplasia and 
systemic lupus. 18 Inflammatory-stimuli-dependent gene, HE- 
RC5, is up-regulated by tumor necrosis factor-alpha or lipo- 
polysaccharide stimulation in late inflammation phase. 19 Several 
studies reported that cigarette smoking adversely affects the 
immune system, particularly lymphocyte development, and 
strongly correlates with inflammation markers. 20,21 However, 
smoking cessation restores the immune levels to normal and 
lowers airway inflammation, and the degree of improvement 
is better when using the anti-inflammatory treatment. 20,22 Ex- 
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Figure 2. Quantitative analysis of the 
mRNA expression levels relative to be- 
ta-actin between the success group and 
the failure group (mean±SD). Graphs 
A-D show the gene expression levels of 
A: herpa cell adhesion molecule 2, B: 
Tropomodulin 1 , C: Ubiquitin specific pe- 
ptidase 18, D: Interferon-induced protein 
44. "indicates statistical significance (p< 
0.05). 



trapolating from our results, smoking cessation could decrease 
the inflammation, protect humans from the development of 
infections and, overall, may strengthen immunity. Moreover, 
changes in the expression of immune-related genes might par- 
ticipate in the improvement of immunity. 

The second group was cytoskeleton-related genes, which 
include TMOD1, EPB49, ANK1 and EPB41. All these genes 
were more highly expressed in the success group. Cigarette 
smoking has been linked to impairment of cytoskeletal integ- 
rity and endothelial organization. 23,24 EPB49 is an actin-bind- 
ing and bundling protein of the eryfhroid membrane skeleton, 
and is abundantly expressed in the human brain, heart, skeletal 
muscle and the lung, 25 thus, it supports the erythrocyte mem- 
brane. Also, ANK1 is a major protein of the erythrocyte mem- 
brane skeleton that forms the link between the membrane 
and the underlying spectrin cytoskeleton. 26 In addition, EBP41, 
as an essential component of the spectrin-based membrane 
skeleton, 27 it plays a key role in determining the morphologi- 
cal and mechanical properties of red cells. Consequently, it is 
reasonable to speculate that improved the membrane regulates 
physical properties of the mechanical stability, and deform- 
ability after smoking cessation, especially in red blood cells. 

Finally, cell cycle and growth-related genes were classified 
as the third group, which contained MAPK3, CCNG2, ANXA3 
and IFI44. MAPK3, CCNG2 and ANXA3 were up-regulated 
and IFI44 was down-regulated in the success group. MAPK3 
induces cell proliferation or differentiation, such as in the 
neuronal regeneration, after injury and acts in whole-body 



development. 28 CCNG2 might promote cell cycle progres- 
sion, and is significantly down-regulated in oral cancer epi- 
fhelia. 29 A recent study reported ANXA3 acts as a positive reg- 
ulator on hepatocyte growth 30 and decreased IFI44 expression 
levels probably allow normal cell growth. 31 Previous cell stud- 
ies have concluded that cigarette smoking has negative effects 
on cytotrophoblast progression through the cell cycle and 
prevents cell growth by inducing significant cytotoxicity. 32,33 
Our results supported these previously described outcomes. 
Furthermore, our findings indicate that, during smoking ces- 
sation, cell proliferation and development will increase due 
to cellular compensation. Although genes involved in medi- 
ating smoking cessation may still be unknown or functional- 
ly uncharacterized, they may also have a role in reduced pre- 
ventable diseases caused by smoking. 

We validated the microarray results obtained with some rep- 
resentative genes considered as a highlight potential biomarkers 
that are participate in different cellular and physiological func- 
tions in a smoking cessation. The 4 genes, i.e., HEPACAM2, 
TMOD1, USP18 and IFI44 were selected because they might 
reduce preventable diseases caused by smoking and manage 
specific physiological activates related to smoking cessation. 
The mRNA expression levels of HEPACAM2 and TMOD were 
up-regulated and USP18 was down-regulated in the success 
group. HEPACAM2 mediates cell- extracellular matrix and 
cell-cell interactions, which are critical for the formation and 
maintenance of cellular architecture, but also for normal bio- 
logical processes, which includes the regulation of cellular 
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adhesion, proliferation, apoptosis, migration and differentia- 
tion. 34 Moreover, its expression is frequently silenced in a va- 
riety of tumor types, because disruption of normal cell adhe- 
sion in the transformed cells may contribute to the enhanced 
migration and proliferation of tumor cells, leading to inva- 
sion and metastasis. 34,35 In addition, TMOD1 is a highly con- 
served capping protein of the pointed ends of the sarcomeric 
actin filaments and erythrocyte membrane, and is involved 
in the architecture of the sarcomere in muscle cells and the 
membrane skeleton in the non-muscle cells. 36 TMODl-null 
mice exhibit a mild anemia with features that resemble he- 
reditary spherocytic elliptocytosis, failure of myofibril assem- 
bly and aborted cardiac development. 37,38 The current results 
support the suggestion that USP18 might be linked to the in- 
nate immune response to viruses and USP18 gene down-reg- 
ulation, which results in an increased antiviral activity against 
many kinds of viruses. 39 Interestingly, USP18 probably has on- 
cogenic properties, and its inhibition may have possibilities for 
cancer therapeutics. 40 Cohort studies have indicated that pre- 
vention and cessation of smoking are the most important de- 
terminants of lung, head and neck cancer. 41,42 Hence, it can be 
assumed that the effects of the increased expression of the 
HEPACAM2 and TMOD1 genes, and the lowered USP18 
gene expression resulting from smoking cessation affects the 
red blood cell membrane and antiviral ability This may also be 
associated with reduced risk of cigarette smoking-related dis- 
eases, especially cancers. 

To the best of our knowledge, this is the first study about the 
effectiveness of smoking cessation on gene expression of smok- 
ers, while further studies are required to delineate the biologi- 
cal response involved in intervention of impairment caused by 
cigarette smoking and development of systemic diseases. This 
microarray analysis demonstrates that smoking cessation of 
heavy smokers represented the changes in the expression of 
at least 299 genes in the circulating blood, and indicates that 
smoking cessation alters the gene expression that is related to 
immunity cytoskeleton, cell adhesion, cell proliferation and 
cancer. The present data provides insights into smoking ces- 
sation-related pathological processes by determining its in- 
fluences at the gene expression level. Gene profiles obtained 
from these results, therefore, could serve as available infor- 
mation regarding the beneficial effects of smoking cessation 
and help elucidate the specific mechanisms that mediate the 
host response to cigarette smoke exposure. 
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